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Methanol  decomposition  on  Pt/Ce02/Zr02  catalyst  is  studied  inside  a  packed  bed  microreactor  in  the 
temperature  range  of  300-380  °C.  The  microreactor  is  fabricated  using  low-temperature  co-fired 
ceramic  (LTCC)  technology,  which  is  well  suited  for  the  production  of  relatively  complex  three- 
dimensional  structures.  It  is  packed  with  2  wt%  Pt— CeC>2  catalyst,  which  is  deposited  onto  Zr02  spher¬ 
ical  particles.  A  ID  mathematical  model,  which  incorporates  diffusion,  convection  and  mass  transfer 
through  the  boundary  layer  to  the  catalyst  particles,  as  well  as  a  3D  computational  fluid  dynamics  model, 
are  developed  to  describe  the  methanol  decomposition  process  inside  the  packed  bed.  The  microreactor 
exhibits  reliable  operation  and  no  catalyst  deactivation  was  observed  during  three  months  of  experi¬ 
mentation.  A  comparison  between  the  ID  mathematical  model  and  the  3D  model,  considering  the  full  3D 
geometry  of  the  microreactor  is  made  and  the  differences  between  the  models  are  identified  and 
evaluated. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Hydrogen  has  been  recognized  as  one  of  the  more  important 
energy  carriers  of  the  future  [1],  Its  suitability  for  miniaturized 
portable  fuel  cell  power  sources  is  mainly  because  of  its  high  energy 
density  on  a  mass  basis  and  microchemical  systems  are  especially 
suitable  for  portable  fuel  cells  due  to  efficient  heat  and  mass 
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transfer  characteristics  and  compact  size  [2].  The  use  of  fuel  cells  for 
generating  electricity  is  showing  promise  since  they  are  an  envi¬ 
ronmentally  friendly,  clean  and  efficient  energy  source  [3]. 

It  has  been  demonstrated  that  microreactors  offer  superior 
operation  as  opposed  to  conventional  fixed-bed  catalytic  reactors 
in  terms  of  methanol  conversion  and  hydrogen  selectivity,  most 
likely  due  to  the  advantage  of  being  able  to  achieve  optimal  process 
conditions  for  the  desired  chemical  transformation,  and  by  taking 
advantage  of  the  large  surface  area  [4].  The  methanol  decomposi¬ 
tion  reaction  is  endothermic  and  is  favored  at  high  temperatures 
[5],  which  necessitates  efficient  heat  transfer  characteristics  of  the 


A  Pohar  et  al.  /  Journal  of  Power  Sources  256  (2014)  80-87 


81 


system.  For  microreactor  fabrication,  low  temperature  co-fired 
ceramic  (LTCC)  technology  is  a  promising  technique  for  producing 
relatively  complex  three-dimensional  multilayer  structures,  by  the 
bonding  of  several  ceramic  layers  [6].  The  integration  of  sensors 
and  actuators  with  microelectromechanical  systems  (MEMS) 
technology  has  been  reported  [7], 

For  hydrogen  storage,  methanol  is  the  most  attractive  hydrogen 
fuel  carrier  due  to  its  high  hydrogen  to  carbon  ratio.  With  the 
exception  of  formic  acid,  methanol  is  the  easiest  hydrocarbon  to  re¬ 
form  and  liquid  methanol  can  readily  be  produced  from  biomass, 
unlike  gasoline  or  diesel  fuel  [8].  It  is  also  much  safer  and  more 
convenient  to  distribute  and  store  hydrogen  chemically  bonded  in 
liquid  form  than  as  a  pure  compound  for  long-distance  transportation 
or  through  gas  pipelines  [9].  Moreover,  methanol  is  sulfur- free  and  the 
absence  of  carbon-carbon  bonds  reduces  the  risk  of  coking  [1  ]. 

Among  the  processes  to  recover  hydrogen  from  methanol  the 
most  frequently  used  is  methanol  steam  reforming,  which  gives  the 
highest  concentration  of  hydrogen  and  the  lowest  concentration  of 
carbon  monoxide  in  the  reformate  stream.  However,  it  requires  an 
overstoichiometric  concentration  of  water  in  the  feed,  which  may 
causes  degradation  of  phosphoric  acid  doped  polybenzimidazole 
high-temperature  PEM  membrane  in  the  fuel  cell.  In  order  to  avoid 
this  undesired  influence,  one  can  use  a  methanol  decomposition 
based  fuel  processor  [10]. 

Pt/Ce02  catalyst  has  been  reported  to  be  among  the  best  for  the 
methanol  decomposition  reaction  [11].  The  addition  of  CeC^  sta¬ 
bilizes  catalyst  operation  by  preventing  carbon  decomposition  on 
the  surface,  which  leads  to  enhancements  in  hydrogen  selectivity 
and  yield  [12].  Pt  and  Pd  based  catalysts  promoted  with  ceria  have 
been  found  to  be  significantly  more  active  than  catalysts,  which  are 
not  promoted.  Furthermore,  the  presence  of  ceria  decreases  the 
temperature  required  for  complete  conversion,  which  was 
explained  as  being  most  likely  due  to  the  ability  of  ceria  to  stabilize 
small  Pt  crystallites  [8]. 

In  this  work,  methanol  decomposition  inside  a  LTCC  micro¬ 
reactor  packed  with  Pt/Ce02/Zr02  catalyst  was  studied.  The  kinetic 
parameters  were  evaluated  and  a  comparison  was  made  between  a 
simplified  ID  model  and  a  3D  model,  which  took  into  account  the 
full  3D  geometry  of  the  micro  reactor. 

2.  Materials  and  methods 

2.1.  Catalyst  preparation 

Mesoporous  ceria  support  having  a  BET  specific  surface  area  of 
132  m2  g-1  was  synthesized  by  using  a  hard  template  method  with 
KIT-6  silica  [13].  The  2  wt%  Pt-Ce02  catalyst  precursor  was  pre¬ 
pared  by  wet  impregnation  of  CeCh  with  a  solution  of  dihydrogen 
hexachloroplatinate  hydrate  (>99.9%,  Sigma- Aldrich,  St.  Louis, 
USA)  in  ethanol.  The  resulting  material  was  dried  overnight  at  60  °C 
in  ambient  air.  A  suspension  of  the  prepared  catalyst  precursor  was 
further  deposited  with  a  method  of  vacuum  evaporation  onto 
350  pm  Zr02  spherical  beads.  The  Zr02  beads,  coated  with  the 
catalyst  precursor,  were  then  dried  at  125  °C  for  2  h  and  subse¬ 
quently  calcinated  at  550  °C  for  3  h  in  a  flow  of  air  to  obtain  a  core¬ 
shell  catalyst.  Before  catalyst  testing,  it  was  reduced  in  situ  at  220  °C 
for  4  h  in  a  flow  of  5  vol%  H2  in  Ar.  Scanning  electron  microscope 
(SEM)  images  of  the  catalyst  were  made.  The  catalyst  BET  surface 
area  was  measured  with  a  MicroMeritics  ASAP2010  instrument. 

2.2.  Microreactor  design 

The  microreactor  was  fabricated  with  LTCC  technology,  as 
described  elsewhere  [14].  Single  ceramic  layers  were  manufac¬ 
tured,  laminated  together  and  fired  in  one  step  at  a  low  firing 


temperature  of  approximately  850  °C.  The  micro  reactor  consisted 
of  inflow  and  outflow  ports  for  fluidic  connections  (</>  2.6  mm),  4 
mounting  ports  (</>  4.5  mm),  and  a  service  port  for  the  insertion  of 
the  catalytic  spheres  (Fig.  la).  A  grate  at  the  entrance  and  exit  of  the 
catalytic  bed  with  76  ducts,  each  with  a  dimension  of  1  x  0.15  mm, 
was  fabricated  for  the  entrapment  of  the  catalytic  particles  and  for 
allowing  the  passage  of  gas.  The  catalytic  beads  were  poured  into 
the  microreactor  through  the  service  hole,  which  was  finally  filled 
with  fiberglass  and  closed  with  graphite  sealing  tape  and  Swagelok 
connection  ports,  which  were  fitted  at  the  inlet  and  outlet  (Fig.  lb). 
The  use  of  spherical  catalyst  particles  for  the  packed  bed  is  ad¬ 
vantageous,  since  closed  areas  inaccessible  to  the  flow  and  stagnant 
zones  are  less  likely  to  exist. 

The  catalytic  bed  had  the  dimensions  33.8  x  3  x  41.4  mm 
(width  x  depth  x  length),  which  yielded  approximately  4.2  cm3 
reaction  chamber  space.  The  total  microreactor  size  was 
63.8  x  6  x  41.8  x  mm. 

2.3.  Process  set-up 

The  process  set-up  is  illustrated  on  Fig.  2.  Laboratory  syringe 
pumps  (TSE  Systems,  Chesterfield,  USA)  with  5  ml  glass  syringes 
were  used  for  methanol  dosage  (of  analytical  grade;  Sigma- 
Aldrich,  St.  Louis,  USA).  An  evaporator,  consisting  of  a  1/16"  stain¬ 
less  steel  coil  capillary  with  an  integrated  heater,  was  used  for  the 
production  of  vapor  methanol.  The  evaporator  was  kept  at  a  con¬ 
stant  temperature  between  120  and  150  °C,  depending  on  the  flow 
rate.  Vaporized  methanol  was  subsequently  fed  into  the  micro¬ 
reactor,  which  was  positioned  inside  an  electric  heater  with  a 
quartz  sand  bath,  thermostated  at  the  operational  temperature. 
Due  to  the  endothermic  reaction,  sufficient  heating  was  provided 
and  the  preset  microreactor  temperature  was  kept  constant  within 
±1°  in  the  temperature  range  of  300-380  °C.  The  gas  from  the 
microreactor  outlet  was  directly  connected  to  a  gas  chromatograph 
for  analysis.  Nitrogen  gas  was  preheated  and  added  to  the  gaseous 
reaction  stream  through  a  mass  flow  controller  and  was  used  as  the 
internal  standard.  Experimentation  was  conducted  at  liquid 
methanol  flow  rates  of  2-70  pL  min-1,  which  corresponds  to  vol¬ 
ume  hourly  space  velocities  of  400-13,700  h-1.  The  volume  hourly 
space  velocity  is  defined  as  the  total  volumetric  flow  rate  to  the 
microreactor  at  standard  temperature  and  pressure,  divided  by  the 
volume  of  the  catalyst  in  the  microreactor. 

2.4.  Gas  chromatographic  analysis 

Vapor  composition  from  the  microreactor  exit  was  analyzed 
with  an  Agilent  7890A  gas  chromatograph  (Agilent  Technologies, 
Santa  Clara,  USA).  Two  capillary  columns:  HP-PLOT-Q 
(30  m  x  0.53  mm  x  40  pm)  and  HP-PLOT-5A 
(30  m  x  0.53  mm  x  50  pm)  with  a  TCD  detector  at  200  °C,  equipped 
with  molecular  sieves,  were  used,  with  helium  as  the  carrier  gas. 
Oven  temperature  was  kept  at  140  °C  for  12  min  after  which  the 
temperature  was  decreased  at  the  rate  of  20  °C  min-1100  °C.  The 
retention  times  of  methanol,  hydrogen,  nitrogen  and  carbon 
monoxide  were  identified  at  5.1,  9.2,  15.1,  15.8,  and  17.1  min, 
respectively.  Quantitative  analysis  of  the  reactants  and  products 
was  enabled  after  the  calibration  with  standard  gas  mixtures. 

3.  Modeling 

3.1.  One -dimensional  model 

In  the  packed  bed  microreactor,  methanol  is  decomposed  to 
carbon  monoxide  and  hydrogen  (CH3OH  -►  CO  +  2H2).  The 
decomposition  reaction  is  endothermic  with  AH298  =  90.7  kj  mol-1 
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Fig.  1.  (a)  The  bottom  layer  of  the  LTCC  microreactor  displaying  the  catalytic  bed  with  a  grate  at  the  inlet  and  outlet,  (b)  The  assembled  LTCC  microreactor  with  metal  fittings,  (c)  The 
model  of  the  microreactor;  units  in  meters,  (d)  A  close-up  of  the  grate  connecting  the  inlet  channel  with  the  reaction  chamber. 


[8].  The  adsorption  equilibrium  of  the  reacting  species  can  be 
written  [11]: 

CH3OH  +  2Pt  <-►  CH30— Pt  +  H-Pt 
H2  +  2Pt  <-►  2H— Pt 
CO  +  Pt  <->  CO-Pt 


The  equilibrium  constants  of  methanol  (Km),  hydrogen  (K h)  and 
carbon  monoxide  (JCc)  are  therefore: 


/Cm  = 


uCH3OuH 

PMeOH^v 


Kh  = 


PhX 


Kc  = 


°_co 

P CO^v 


(1) 


Fig.  2.  The  experimental  set-up;  (a)  mass  flow  controllers  (MFCs),  (b)  heater  for  the 
microreactor,  (c)  the  control  unit  for  MFCs,  (d)  the  control  unit  for  additional  heaters 
(evaporator,  nitrogen  inlet,  inlet  to  the  gas  chromatograph),  (e)  laboratory  syringe 
pump,  (f)  evaporator,  (g)  nitrogen  inlet,  (h)  thermostated  microreactor  inside  the 
heater,  (i)  gas  chromatograph,  (j)  temperature  indicator,  (k)  computer  analysis 
software. 


where  ^ch3o^h^co  are  the  concentrations  of  the  species  bound  to 
the  catalyst  and  dv  is  the  concentration  of  the  vacant  sites.  P;  is  the 
partial  pressure  of  each  species. 

The  dominant  dissociation  pathway  of  methanol  involves  the 
cleavage  of  the  0— H  bond  [15],  while  the  abstraction  of  the 
hydrogen  atom  from  the  methyl  group  was  found  to  be  the  rate 
determining  step  of  methanol  decomposition  [11,16]: 

CH30— Pt  +  H-Pt  H2CO— Pt  +  H2  +  Pt  k 


where  k  is  the  kinetic  constant.  The  derived  rate  equation  of  the 
limiting  step,  considering  0Ch3o  +  #h  +  #co  +  0V  =  1,  is: 

r  =  _ /cKmPMeOH _ 

(1  +  XmX“1/2PMeOHPH21/2  +  K-'/2PHy2  +  KcPco) 2 

Due  to  the  small  inner  reactor  dimensions  and  short  diffusion 
pathways,  the  operation  of  the  miniaturized  reactor  was  described 
with  a  one-dimensional  isothermal  model  with  x  as  the  length 
coordinate.  Diffusion,  convection  and  mass  transfer  through  the 
boundary  layer  to  the  catalyst  particles  were  included  in  the  mass 
balance  of  bulk  flow.  For  low  values  of  the  Reynolds  number 
{Re  <  1 ),  the  dispersion  coefficient  is  equal  to  the  molecular  diffu¬ 
sion  coefficient  [17]. 


dc^  d(t'Cj)  _  (1  -e) 

dx  £ 


dx  \  !  dx  ' 


kc,ia(ci  Ccat,i)  —  0 


(3) 


Index  i  represents  methanol,  nitrogen,  hydrogen  and  carbon 
monoxide.  The  Danckwerts  boundary  conditions  are  as  follows: 


_DXi 

U'dx 


x=0+ 


vCi\o-  ~  vCi\o+  x  =  0 


dc,/dx  =  0  x->oo 

where  Dj  is  the  diffusion  coefficient  of  each  species,  q  is  the  bulk 
concentration,  ccat,i  is  the  concentration  at  the  catalyst  surface,  v  is 
the  interstitial  velocity,  and  a  is  the  surface  area  per  void  volume. 
The  mass  balance  at  the  catalyst  surface  is: 
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kc.i+cat.i  -  Cj)  =  W  (4) 

where  v\  is  the  coefficient  of  stoichiometry  and  /<c,i  is  the  mass 
transfer  coefficient.  Vapor  viscosities  of  the  non-polar  gases 
involved  in  the  reaction  were  predicted  with  the  Stiel  and  Thodos 
correlation  [18]  from  their  critical  temperatures  (TCI)  and  pressures 

(Pc,*): 


Cj 

E  cf 


(11) 


Similarly,  the  gas  mixture  density  was  calculated  considering 
ideal  gas  behavior: 


Pmix  ~  MjXj 


(12) 


Vi  =  4.6  x  10-4jVr 


My2pc2f 

rVe’ 

c,i 


where  P  is  the  pressure,  R  is  the  gas  constant,  and  T  is  the  tem- 
(5)  perature.  The  pressure  drop  for  spherical  packing  in  a  packed  bed 
reactor  was  calculated  with  the  Ergun  equation: 


where  Mj  is  the  molar  mass  in  g  mol  1  and 
Ni  =  3.4  x  10-4Tr0-94  Tr  <  1.5 


Nj  =  1.778  x  (4.58Tr  j  -  1.67)°  625  Tr  >  1.5 


Tr,i  =  T/TcJ  (6) 

Viscosities  of  polar  gases  were  calculated  from  the  method  of 
Chung  et  al.  [19  ,  which  accounts  for  the  molecular  shape  and  the 
polarity  of  polar  gases: 


Vi 


40.785 


V2C fovj 


(7) 


in  which  case  the  viscosity  is  in  pP,  the  molar  mass  is  in  g  mol-1,  Vc>1 
is  the  critical  volume  in  cm3  mol-1,  and  Qv,i  is  the  viscosity  collision 
integral,  defined  as  [20]: 


/  -0.14874  / 

QVji  =  1.16145 +0.52487  exp (-0.773207^ 
+  2.16178  exp(-2.437877/[) 


T(  =  1. 25937V  i  (8) 

Eq.  8  is  valid  in  the  range  0.3  <  T*  <  100,  with  an  average  de¬ 
viation  of  only  0.064%  [21  ].  The  factor  FCiI-  accounts  for  the  molecular 
shapes  and  polarities  of  dilute  gases: 

FC  j  =  1  -  0.2756w,  +  0.059035+  +  K; 


dP 

dx 


G  /I  -  e\  /150(1  -e)v 

PDp  17“ j  V  D~P 


1.75G 


(13) 


which  takes  into  account  the  characteristics  of  the  packing  (void 
fraction  e ,  particle  diameter  Dp),  the  characteristics  of  fluid  prop¬ 
erties  (viscosity  v ,  density  p),  and  the  characteristics  of  the  oper¬ 
ating  parameters  (mass  velocity  G).  The  pressure  was  calculated 
from  the  outlet  of  the  microchannel  (at  1  atm)  to  the  inlet  and 
iterated  along  with  the  equations  of  state.  The  packing  character¬ 
istics  were  considered  invariant  with  respect  to  reactor  length  and 
the  operation  was  assumed  to  be  steady-state  and  isothermal. 

The  mass  transfer  coefficients  (/+/)  were  calculated  from  the 
correlation  for  the  Sherwood  numbers  (Shi)  through  an  expression 
proposed  by  Comiti  et  al.  [22]  at  each  longitudinal  position  of  the 
channel: 


Shi 


3.66  +  0.101  ^1  +Er(T_  ^ 
x  (l+^Sc,1/3exp(-ll^AYXe+)))  xXe+48Sc+ 


kri  = 


ShjDi 

“DT 


(14) 


The  equation  is  valid  for  laminar  and  turbulent  flow  at 
0.03  <  Re  <  105.  t  represents  tortuosity,  which  is  the  ratio  of  the 
actual  length  of  the  flow  path  to  the  straight-line  distance  between 
the  ends  of  the  flow  path,  and  is  therefore  tc/2  for  spherical  parti¬ 
cles.  Y  is  the  particle  aspect  ratio,  defined  in  the  direction  normal  to 
the  flow,  and  is  equal  to  1  for  spheres.  Sq  is  the  dimensionless 
Schmidt  number  (Scz  =  p/pD/)  and  Xew  is  the  dimensionless  wall 
energetic  criterion,  which  can  be  calculated  for  porous  media: 


Ki  =  0.0682  +  4.704 N0H,i/Mi 


Mr,i 


=  131.3 


Mi 

(Vc,iTCli) 


1/2 


(9) 


where  <+■  is  the  acentric  factor,  /q  is  a  correction  for  highly  polar 
substances,  jUr.i  is  a  dimensionless  dipole  moment,  ^  is  the  dipole 
moment  in  debyes,  and  N0H,i  is  the  number  of  -OH  groups. 

The  viscosity  of  the  gaseous  mixture  is  a  function  of  pressure 
and  composition  and  varies  along  the  microreactor  length.  It  was 
calculated  via  the  following  expression,  which  is  accurate  at  low 
pressures,  at  each  segment  of  the  microreactor: 


Ex,-4,-(M,-)1/2 


(10) 


where  x2-  is  the  mole  fraction  of  component  i  with  viscosity  Vi  and 
molar  mass  Mi  and  is  calculated  from: 


Xew  =  64Re2  (l  +  1.21  x  10  ~2Re)  (15) 

The  Reynolds  number,  for  fluid  flow  in  a  packed  bed  of 
approximately  spherical  particles,  is: 


Re 


pVsDp 

V(\-e) 


(16) 


where  Vs  is  the  superficial  velocity  of  the  flow. 

The  diffusion  coefficients  in  multicomponent  gas  mixtures  of 
the  species  under  consideration  were  calculated  from  the  binary 
diffusion  coefficients  through  the  following  relation  [23]: 


(17) 


where  Da  is  the  effective  diffusion  coefficient  of  gas  A  with  respect 
to  the  total  gas  mixture,  while  DAb,  Dac,  etc.  are  the  binary  diffusion 
coefficients. 
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The  diffusion  coefficients  of  binary  gas  mixtures  were  calculated 
with  a  diffusion  volume  correlation  method,  developed  by  Fuller  j * 

et  al.  [24]: 


-  (  P^i  ^2  Dikdi 


(25) 


dab 


IQ-8!1  75(1  /Ma  +  1  /Mb)1/2 


P 


(18) 


in  which  case  Dab  is  in  the  units  of  cm2  s_1,  P  is  the  pressure  in  atm, 
and  Vj  is  a  special  atomic  diffusion  volume  parameter,  provided  in 
the  reference  for  each  constitutive  atom  of  the  gas. 

The  velocity  increase  of  the  gaseous  mixture  due  to  the  reaction 
was  calculated  considering  ideal  gas  behavior.  Mass  fractions  (w/) 
and  partial  pressures  (P,-)  were  calculated  from  molar  fractions: 


dk  =  Vxfc  +  l((xfc-wfc)VP)  (26) 

where  Dik  are  the  multicomponent  Fick  diffusivities  and  dk  is  the 
diffusional  driving  force  acting  on  species  k.  The  kinetic  expression 
is  presented  in  Eq.  (2).  Mass  transfer  limitations  were  not  included 
due  to  the  information  obtained  from  the  results  of  the  1 D  model. 
The  boundary  conditions  are  the  Danckwerts  boundary  conditions: 
Inlet: 

ji\y=0+  —  ^inPwily=0-  —  vinP(x)i\y=0+  (27) 


Wf 


Vin  is  the  inlet  normal  velocity  calculated  from  the  flow  rate  and  the 
(19)  inlet  area  (An). 

Outlet: 


Pi  =  *iP 


(20)  j.  =  0;  P  =  1  atm 


(28) 


The  vapor  properties  were  evaluated  at  every  position  in  the 
microreactor  and  were  corrected  with  the  pressure  field  in  the 
iterative  solving  process.  Finite  differences  were  used  for  the  nu¬ 
merical  discretization  with  the  inclusion  of  a  pseudo  time  step  for 
the  treatment  of  the  second  derivative  in  the  mass  balance  equa¬ 
tions.  The  discretization  of  the  computational  domain  was  refined 
until  no  grid  dependence  was  observed.  A  60  x  30,000 
(length  x  time)  discretization  was  used  with  the  initial  conditions 
the  same  as  the  inlet  boundary  conditions.  The  solution  was  iter¬ 
ated  until  it  converged  to  steady  state  with  a  desired  accuracy  level. 
Matlab  was  used  for  performing  the  calculations. 


The  no  flux  and  no-slip  boundary  conditions  were  used  on  the 
remaining  boundaries.  The  permeability  of  the  packed  bed  (k)  was 
calculated  with  the  correlation  published  by  Rumpf  and  Gupte  [25]: 


£5’5  - 

=  5^ 


(29) 


Several  grid  sizes  were  tested  to  ensure  that  the  results  were 
grid  insensitive  and  573,511  tetrahedral  elements  were  finally  used. 
The  convergence  criteria  for  the  normalized  residuals  for  each 
variable  were  restricted  to  less  than  10~6. 


3.2.  Three-dimensional  CFD  model 

Additionally,  a  3D  finite  element  steady-state  model  was  devel¬ 
oped.  The  microchannel  geometry  was  modeled,  meshed,  and  solved 
with  COMSOL  Multiphysics.  The  model  consisted  of  an  inlet  and 
outlet  channel,  and  a  catalytic  reaction  chamber  (packed  bed  with 
porous  media  flow),  where  the  reaction  took  place.  The  pressure  and 
the  velocity  were  continuous  over  the  interface  between  the  free  and 
porous  domain.  Gravitational  effects  were  neglected.  Continuity  and 
momentum  balance  equations  (Navier- Stokes  equations)  were  used 
for  modeling  the  compressible  laminar  flow: 

V(pu)  =  0  (21) 

p(  U-V)u  =  -VP  +  V-(u(vu  +  (Vu)T))  +  (_^v-u  j  (22) 

where  u  is  the  velocity  vector.  The  Brinkman  equation  was  used  for 
modeling  flow  through  a  porous  matrix: 

'((“V>")  -  -VP  +  V- (I  (,u(vu  +  (Vu)T))  +  (-fvu)) 


(23) 

where  kp  is  the  permeability  of  the  porous  medium. 

The  Maxwell— Stefan  diffusion  model  was  used  for  the 
description  of  mass  transfer: 

V-ji  +  p(u-V)Wi  =  rf  (24) 


4.  Results  and  discussion 

The  reaction  chamber  was  filled  with  approximately  1.16  x  105 
Pt/Ce02/Zr02  catalyst  spherical  particles  with  an  average  diameter 
(Dp)  of  350  pm  and  a  total  mass  of  12.5  g.  The  total  amount  of  the 
catalyst  on  the  particles,  obtained  by  weighting  coated  and  un¬ 
coated  catalyst  particles,  is  83  mg.  The  calculated  void  volume  (e)  of 
the  chamber  is  0.38  and  the  permeability  (k)  is  1.07  x  10-10  m2 
(from  Eq.  (29)).  BET  analysis  determined  a  specific  surface  area  of 
0.2  m2  g-1.  Fig.  3  presents  the  SEM  image  of  the  coated  core-shell 
catalyst.  The  average  thickness  of  the  shell  is  approximately  5  pm. 
Energy-dispersive  X-ray  spectroscopy  (EDS)  performed  on  the 
catalyst  particles  identified  38.7  wt%  O,  34.4  wt%  Zr,  25.1  wt%  Ce, 
and  1.73  wt%  Pt. 

Numerical  simulations  show  that  the  rate  expression  (Eq.  (2))  in 
steady  state  can  be  reduced  to: 

r  m  kDI(MeKCQ PMe^co  (30) 

The  natural  logarithms  of  kDI<MeI<^  values,  obtained  with 
regression  analysis  at  each  experimental  temperature  (300,  317, 
340,  355,  380  °C),  in  dependency  to  reciprocal  temperature,  give  a 
linear  correlation  (Fig.  4,  Table  1). 

The  results  at  380  °C  and  a  flow  rate  of  10  pL  min-1  liquid 
methanol  with  10  wt%  N2  show  an  increase  of  the  mean  velocity  in 
the  cross-section  of  the  reaction  chamber  (yz  plane)  from  0.095  to 
0.11  m  s-1  (15.8%)  due  to  the  change  in  gas  density,  which  is  a 
consequence  of  the  reaction  stoichiometry,  producing  3  mol  of  the 
products  from  1  mol  of  the  reactant.  The  pressure  drop  inside  the 
microreactor  was  mostly  linear  with  a  pressure  difference  of  125  Pa. 
The  velocity  distribution  across  the  chamber  cross-section  is 
mainly  uniform,  apart  from  some  variation  at  the  entrance  and  exit 
of  the  chamber,  due  to  the  bottom  positioned  inlet  and  outlet  grate. 
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Fig.  3.  SEM  image  of  the  Zr02  spherical  beads  (core)  with  2  wt%  Pt/Ce02  catalyst  shell. 
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Fig.  4.  The  natural  logarithm  of  the  kinetic  parameters  as  a  function  of  reciprocal 
temperature. 


Table  1 

The  values  of  the  kinetic  expression,  which  consists  of  the  kinetic 
and  equilibrium  constants,  obtained  by  regression  analysis. 


T(°C) 

kDI<MeI<co  (mol  s  1  m  3  Pa) 

300 

4.4  x  102 

317 

1.5  x  103 

338 

8.7  x  103 

355 

2.9  x  104 

380 

2.2  x  105 

Fig.  5  presents  the  velocity  magnitude  on  3  yz  and  1  xy  planes.  The 
xy  plane  is  positioned  0.2  mm  above  the  inlet  and  outlet  grate.  The 
color  range  is  set  to  0-0.02  m  s-1  for  a  better  view  of  the  velocity 
magnitude  inside  the  reaction  chamber.  On  the  first  yz  plane  at  the 
top  of  the  chamber  a  lower  velocity  at  the  entrance  can  be  seen;  this 
is  similar  at  the  exit  (not  shown  for  a  better  view  of  the  exit  grate). 
The  highest  velocity  1.237  m  s-1  is  located  inside  the  outlet  grate. 
The  red  lines  present  the  fluid  flow  streamlines,  which  portray  an 
evenly  distributed  flow. 

Fig.  6  presents  the  comparison  between  the  velocity  magnitudes 
of  the  ID  model  and  the  3D  model.  For  the  3D  model,  the  velocity 
magnitude  on  3  lines  at  the  middle  of  the  reaction  chamber, 
running  along  the  length  of  the  reactor,  at  different  depths  were 
chosen  (3D  centerline  -  through  the  middle  of  the  depth,  3D  bot¬ 
tom  1/4  -  at  the  bottom  quarter  of  the  chamber,  and  3D  top  1/4  -  at 
the  top  quarter  of  the  chamber).  The  higher  velocities  at  the 
entrance  and  exit  can  be  observed  in  the  bottom  1/4  of  the  chamber 
due  to  the  bottom  positioned  inlet  and  outlet  grates  with  a  smaller 
cross  section  area.  In  contrast,  the  velocity  at  the  top  1/4  of  the 
chamber,  which  is  located  on  the  other  side  of  the  inlet  and  outlet 
grate,  is  lower  at  the  entrance  and  exit  of  the  reaction  chamber. 

Fig.  7  presents  the  methanol  mass  fraction  distribution.  A  sur¬ 
face  plot  is  presented  due  to  the  small  variance  of  methanol  mass 
fraction  in  the  direction  of  depth.  The  inlet  mass  fraction  is  0.9 
methanol  and  0.1  N2,  which  was  used  as  the  standard  for  the  gas 
chromatography  analysis.  The  catalyst  particles  are  only  present 
inside  the  reaction  chamber,  whereas  a  substantial  lowering  of  the 
methanol  mass  fraction  (and  an  increase  of  CO  and  H2)  can  be 
observed  inside  the  inlet  channel,  due  to  the  fast  diffusion  process, 
which  works  to  level  out  the  concentration  difference.  The  mass 
fraction  of  methanol  entering  the  reaction  chamber  is  already 
reduced  to  around  0.4.  The  mass  fraction  graphs  of  CO  and  FI2  are 
presented  in  Supplementary. 

The  mass  fractions  of  the  species  constituting  the  gas  inside  the 
reaction  chamber  are  presented  on  Fig.  8.  The  3D  model  values  are 
the  calculated  average  values  across  the  chamber  cross-section  (yz 
plane).  The  mass  fraction  of  methanol  is  lower  at  the  entrance  to 
the  reaction  chamber  due  to  the  fast  diffusion,  and  is  accurately 
described  by  the  inlet  Danckwerts  boundary  condition.  Similarly, 
the  mass  fractions  of  H2  and  CO  at  the  entrance  are  already  0.06  and 
0.38.  The  3D  results  predict  a  lower  conversion  compared  to  the  ID 
model  which  is  expected,  since  they  describe  the  deviation  from 
ideal,  ID  operation,  due  to  localized  velocity  and  concentration 
variations.  A  20%  increase  of  the  value  of  the  kinetic  expression  is 
required  for  the  3D  model  to  overlay  the  ID  model  solution. 
Modeling  results  from  the  1 D  model  also  revealed  that  there  was  no 
mass  transfer  resistance  due  to  fast  molecular  diffusion  to  the 
catalyst  surface. 

Fig.  9  shows  the  mass  fractions  in  dependence  of  the  residence 
times.  The  results  of  the  model  were  obtained  by  setting  the  flow 
rates  of  liquid  methanol  between  8  and  1500  pL  min-1,  at  380  °C, 
and  with  10  wt%  N2.  Together  with  the  added  nitrogen  gas,  this 
corresponds  to  inlet  gas  mixture  velocities  of  0.2—37.3  mL  min-1 
and  a  maximum  residence  time  of  8  s.  The  residence  time  i  was 
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Fig.  5.  The  velocity  magnitude  inside  the  packed  bed  microreactor.  The  color  bar  is  for  0-0.02  ms1.  The  red  lines  on  the  graph  represent  the  fluid  streamlines.  (For  interpretation 
of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


Fig.  6.  The  velocity  magnitude  along  the  length  of  the  microreactor  for  the  1 D  model 
and  the  3D  model.  The  velocity  magnitude  of  the  3D  model  is  presented  on  lines 
running  through  the  middle  of  the  microreactor  along  the  length  at  3  different  height 
positions. 


calculated  from  the  void  volume  of  the  microreactor  and  the  flow 
rate.  The  experimental  values  were  obtained  with  a  liquid  meth¬ 
anol  flow  rate  of  10,  40,  and  70  pL  min-1.  Gas  constitution  at  the 
lowest  experimental  flow  rate  was  5.10  wt%  methanol,  71.6  wt%  CO, 
9.25  wt%  H2,  and  14.05  wt%  N2.  With  the  acquired  kinetic  constants, 
a  good  agreement  was  found  between  the  modeling  results  and  the 
experimental  data.  At  the  experimental  flow  rate  of  10  pL  min-1, 
94%  methanol  conversion  was  achieved. 

The  catalyst  exhibited  no  deactivation  after  three  months  of 
intensive  experimentation  and  EDS  analysis  did  not  detect  any 
carbon  deposited.  Matolin  et  al.  [26  found  that  with  Pt/ceria  cat¬ 
alysts,  a  complete  removal  of  carbon-containing  intermediates 
occurs  above  500  K.  The  effect  was  explained  by  the  reaction  of  the 
intermediates  with  oxygen,  provided  by  reverse  spillover  from 
ceria. 

5.  Conclusions 

A  packed  bed  microreactor,  fabricated  with  LTCC  technology,  is 
used  for  studying  methanol  decomposition.  With  the  microchannel 
design  the  velocity  distribution  achieved  inside  the  packed  bed  is 
almost  homogenous.  2  wt%  Pt-Ce02  catalyst  is  a  very  selective  and 
active  methanol  decomposition  catalyst.  The  catalyst  exhibits 


▲  0.9 


Fig.  7.  Methanol  mass  fraction  at  380  °C  and  the  liquid  methanol  flow  rate  of  10  pL  min  \  The  geometric  values  are  in  meters. 
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Fig.  8.  Modeling  results  of  the  mass  fractions  of  the  species  along  the  length  of  the 
reaction  chamber  at  380  °C,  a  liquid  methanol  flow  rate  of  10  pL  min-1  and  10  wt%  N2. 
The  3D  results  are  presented  as  average  values  across  the  cross-section  (yz  plane). 


Fig.  9.  The  mass  fractions  obtained  by  the  ID  model,  3D  model  and  the  experimental 
results  at  380  °C  at  various  residence  times. 

excellent  resistance  to  surface  poisoning  by  formation  of  stable 
carbides,  due  to  the  reaction  with  oxygen,  supplied  by  ceria. 
Methanation  side-reaction  does  not  occur. 

An  accurate  solution  is  obtained  by  the  3D  model  and  the  de¬ 
viation  of  the  microreactor  from  ideal  ID  behavior  is  due  to  non¬ 
ideal  flow  and  concentration  distribution.  Modeling  results  show 
a  significant  concentration  change  taking  place  inside  the  inlet 


channels,  even  though  the  reaction  only  takes  place  inside  the 
catalytic  bed.  Optimization  of  the  microreactor  geometry  would  be 
achieved  by  positioning  the  inlet  and  outlet  grates  in  the  flow  di¬ 
rection.  This  way,  the  gaseous  flow  would  not  undergo  as  many 
directional  changes,  which  would  cause  a  smaller  pressure  drop 
and  would  provide  a  uniform  velocity  distribution.  The  operation  of 
the  reformer  would  approach  ideal  behavior,  as  is  described  by  the 
ID  model. 
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